COMMENT ⊗   VALID 00004 PAGES
C REC  PAGE   DESCRIPTION
C00001 00001
C00002 00002	BEGIN "BETERM - COMPUTES BEJCZY COEF. FROM PHYSICAL PARAMETERS"
C00004 00003	⊃ COMPUTED TERMS
C00006 00004	⊃ TYPE OUT THE RESULTS
C00008 ENDMK
C⊗;
BEGIN "BETERM - COMPUTES BEJCZY COEF. FROM PHYSICAL PARAMETERS"

DEFINE ⊃="COMMENT",CRLF="('15&'12)";

⊃ THE FOLLOWING ARRAYS CONTAIN THE PHYSICAL PARAMETERS;

PRELOAD_WITH  513610.0, 259780.0, 258787.0, 101390.0,  41460.0,  33686.0;
OWN REAL ARRAY M[1:6];

PRELOAD_WITH 0.0,-4.28,0.0,0.0,0.0,0.0; 
OWN REAL ARRAY Y[1:6];

PRELOAD_WITH -1.69,-1.67,-18.5,1.323,2.82,-3.09;  
OWN REAL ARRAY Z[1:6];

PRELOAD_WITH  0.0     , .96329@7, .79694@8, 633330.0, 440346.0, 422775.0;
OWN REAL ARRAY IX[1:6];

PRELOAD_WITH  .54944@7, .20459@7, .79694@8, 174440.0, 414080.0, 395629.0;
OWN REAL ARRAY IY[1:6];

PRELOAD_WITH  0.0     , .77522@7,  39167.0, 505560.0,  38837.0,  39335.0;
OWN REAL ARRAY IZ[1:6];

PRELOAD_WITH 0.0,6.050,0.0,0.0,0.0,10.281;
OWN REAL ARRAY R[1:6];

PRELOAD_WITH  1.5@8, 4.0@8, 100000.0, 1.0212@7, 1.0212@7, 848093.0;
OWN REAL ARRAY SHAFT[1:6];

INTEGER DUM;
⊃ COMPUTED TERMS;

REAL B1,B2,B3,B4,B5,B13,B14,B15,B17,B18,B19,B20,D1,D2,D3,A,G;
EXTERNAL REAL PROCEDURE SQRT(REAL X);

G  ← 3.859198*10↑(-4);
B1 ← SHAFT[1] + IY[1] + 2.0*M[2]*Y[2]*R[2] + (M[2]+M[3]+M[4]+M[5])*R[2]↑2
	+ 0.5*(IX[2]+IZ[2]);
B2 ← IY[3];
B3 ← 2.0*(M[3]*Z[3]-M[4]*Y[4]);
B4 ← SHAFT[3]+M[3]+M[4]+M[5]+M[6];
B5 ← 2.0*M[5]*Z[5] + 2.0*M[6]*(Z[6]+R[6]);
B13 ← SHAFT[2] + IY[2] + IX[3];
B14 ← SHAFT[4] + IY[4];
B15 ← IX[5] + IX[6] + M[6]*R[6]*(R[6]+2.0*Z[6]);
B17 ← SHAFT[5] + IY[5];
B18 ← IX[6] + M[6]*R[6]*(R[6]+2.0*Z[6]);
B19 ← SHAFT[6];
B20 ← IZ[6];
D1 ← M[2]*Z[2] + M[3]*Z[3] - M[4]*Y[4];
D2 ← B4-SHAFT[3];
D3 ← 0.5*B5;

⊃ TYPE OUT THE RESULTS;

OPEN(1,"DSK",0,0,2,120,DUM,DUM);
ENTER(1,"BETERM.DAT",DUM);
SETFORMAT(10,5);
OUT(1,CRLF&CRLF&"COMPUTED TERMS:"&CRLF);
OUT(1,"B1 = "&CVE(B1)&CRLF);
OUT(1,"B2 = "&CVE(B2)&CRLF);
OUT(1,"B3 = "&CVE(B3)&CRLF);
OUT(1,"B4 = "&CVE(B4)&CRLF);
OUT(1,"B5 = "&CVE(B5)&CRLF);
OUT(1,"B13 = "&CVE(B13)&CRLF);
OUT(1,"B14 = "&CVE(B14)&CRLF);
OUT(1,"B15 = "&CVE(B15)&CRLF);
A ← B17 + B18;
OUT(1,"B17+B18 = "&CVE(A)&CRLF);
A ← B19 + B20;
OUT(1,"B19 + B20 = "&CVE(A)&CRLF);
OUT(1,"D1 = "&CVE(D1)&CRLF);
OUT(1,"G ="&CVE(G)&CRLF);
OUT(1,"-G ="&CVE(-G)&CRLF);
A ← G*D2;
OUT(1,"G * D2 = "&CVE(A)&CRLF);
A ← -1.0*G*D3;
OUT(1,"-G * D3 = "&CVE(A)&CRLF);
A ← 2.0*(R[6]+Z[6]);
OUT(1,"2(R6+Z6) = "&CVE(A)&CRLF);
A ← R[6]*(R[6]+2.0*Z[6]);
OUT(1,"R6(R6+2Z6) = "&CVE(A)&CRLF);
RELEASE(1);

END